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ABSTRACT 

Hypoxia inducible factor (HIF) up-regulates the tran- 
scription of a few hundred genes required for the 
adaptation to hypoxia. This restricted set of 
targets is in sharp contrast with the widespread dis- 
tribution of the HIF binding motif throughout the 
genome. Here, we investigated the transcriptional 
response of GYS1 and RUVBL2 genes to hypoxia 
to understand the mechanisms that restrict HIF 
activity toward specific genes. GYS1 and RUVBL2 
genes are encoded by opposite DNA strands and 
separated by a short intergenic region ( 1 kb) that 
contains a functional hypoxia response element 
equidistant to both genes. However, hypoxia 
induced the expression of GYS1 gene only. 
Analysis of the transcriptional response of 
chimeric constructs derived from the intergenic 
region revealed an inhibitory sequence whose 
deletion allowed RUVBL2 induction by HIF. 
Enhancer blocking assays, performed in cell 
culture and transgenic zebrafish, confirmed the 
existence of an insulator element within this 
inhibitory region that could explain the differential 
regulation of GYS1 and RUVBL2 by hypoxia. 
Hence, in this model, the selective response to HIF 
is achieved with the aid of insulator elements. This 
is the first report suggesting a role for insulators 
in the regulation of differential gene expression in 
response to environmental signals. 



INTRODUCTION 

A large number of biochemical reactions require oxygen 
as a substrate and metazoa metabolism is largely depend- 
ent on oxidative phosphorylation. At the cellular level, 
the unbalance between oxygen demand and supply 
(hypoxia) results in the activation of a specific gene 
expression program aimed at increasing oxygen delivery 
and reducing its consumption through metabolic 
reprogramming. This transcriptional response is mostly 
mediated by an evolutionarily conserved family of tran- 
scription factors termed hypoxia inducible factors (HIFs), 
that belong to the basic helix-loop-helix superfamily (1). 
HIFs are heterodimers of a constitutive beta subunit 
(HIFp also known as ARNT), that partners with several 
factors and an alpha subunit (HIFa), whose stability (2) 
and transcriptional activity (3) is regulated by oxygen. 
Under hypoxia, HIFa subunits avoid degradation, bind 
to the constitutively expressed beta subunits and the 
heterodimers translocate to the nucleus where they bind 
to the RCGTG motif within the regulatory regions of 
target genes to promote their transcription (4-6). Several 
works have identified individual HIF targets that, taken 
together, account for the metabolic adaptation and induc- 
tion of angiogenesis observed under hypoxia (7). To gain 
insight into the full range of cellular adaptations to 
hypoxia, several groups recently attempted the global 
identification of HIF-targets (5,6,8-11). Interestingly, all 
these works coincide in that only a few hundred, out of all 
the genes containing RCGTG motifs, are regulated by 
hypoxia. Thus, as it is the case for other transcription 
factors (12), HIF binds only to a small proportion of the 
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potential binding sites (5,6,8-11). The basis for this select- 
ivity is incompletely understood, but several mechanisms 
have been proposed. Among them, the cooperation with 
other transcription factors, have been well characterized in 
some instances (13,14). In the case of HIF, requirement of 
functional HNF-4 (15), AP-1 (16), GATA-2 (16) or ETS 
(17,18) sites for proper hypoxic induction of selected 
targets have been described. In agreement with these 
single locus studies, global analysis of HIF binding sites 
by means of experimental (5) and computational methods 
(11) showed the existence of overrepresented transcription 
factor binding sites (TFBS) in close proximity to the 
hypoxia response element (HRE) that might account for 
factors cooperating with HIF. However, the experimental 
characterization of the role of these TFBS in the regula- 
tion of HIF targets by hypoxia is yet to be determined. 
Thus, the cooperation between HIF and other factors 
could contribute to the target selectivity, but it is yet 
unclear to what extent this mechanism explains the 
observed pattern of targets. 

A further mechanism that could dictate the target se- 
lectivity is the accessibility to the TFBS. Histone modifi- 
cations alter the structure of chromatin and hence the 
availability of the underlying nucleotide sequence for the 
binding of transcription factors (19). In addition, DNA 
methylation can preclude the binding of specific transcrip- 
tion factors (20-22). In this regard, a recent study ad- 
dressed the cell-type specificity in response to hypoxia 
and concluded that only those loci that were transcription- 
ally active under basal (normoxic) conditions were permis- 
sive to HIF-regulation (8). However, although these 
results explain most of the intercellular variation in the 
hypoxic transcriptome, it is clear that an additional layer 
of regulation is required, as only a small fraction of all the 
active genes under basal conditions were induced by 
hypoxia in any of the cell lines studied. 

Finally, insulators are included among the regulatory 
mechanisms employed by eukaryotes to ensure specific 
patterns of gene expression and as such, they could be 
involved in the selection of genes to be activated by HIF 
in response to hypoxia. Insulators are defined as DNA 
elements that partition chromatin into independent tran- 
scriptional domains, thereby contributing, in combination 
with additional epigenetic mechanisms, to the tight 
control of gene expression (23) and to the nuclear struc- 
ture and dynamic organization (24). These elements have 
been functionally described, according to their ability to 
block the spread of heterochromatin (barrier function) 
into adjacent loci and to prevent the promiscuous inter- 
action of distal enhancers with proximal promoters, when 
placed in between (enhancer blocking function) (25). The 
role of insulators in the determination of linage-specific 
patterns of gene expression is well characterized, as 
illustrated in the chicken P-globin locus (26). However, 
the participation of these elements in the target discrimin- 
ation by acutely activated transcription factors, in 
response to environmental factors, has not been reported. 

Here, we investigate the mechanism that restricts HIF 
activity toward specific genes by the study of a locus that 
we consider paradigmatic, the GYS1/RUVBL2 genomic 
region from the human genome. The GYS1 gene was 



recently described as a novel hypoxia-inducible gene and 
a functional HRE was identified upstream its promoter 
(27). Interestingly, very close to GYS1, but encoded by 
the opposite DNA strand, is located the RUVBL2 gene. 
In spite of the location of the HRE between both genes, 
we found that only GYS1, but not RUVBL2, was induced 
in response to hypoxia. The lack of RUVBL2 response to 
hypoxia was not due to epigenetic silencing of its promoter 
as it showed a substantial transcriptional activity and the 
level of RUVBL2 mRNA was comparable with that of 
GYS1. Instead, the analysis of different reporter con- 
structs derived from the intergenic GYS1/RUVBL2 
sequence, revealed an inhibitory region, located between 
the HRE and the RUVBL2 minimal promoter, that pre- 
vented the induction of RUVBL2 by HIF. Removal of 
this region allowed the up-regulation of the RUVBL2 
promoter upon HRE activation. This result hinted the 
existence of an enhancer blocking element within the in- 
hibitory region that prevented the interaction between the 
HRE and the RUVBL2 promoter. By means of specific 
enhancer blocking assays (EBA), performed in cell culture 
and using heterologous constructs in transgenic zebrafish, 
we confirmed the existence of an insulator element within 
this locus that could explain the differential regulation of 
GYS1 and RUVBL2 by hypoxia. Altogether, our results 
suggest that HIF selectivity is achieved, at least in this 
locus, by an insulator element that prevents the activity 
of the HRE/HIF complex on the RUVBL2 promoter. 

MATERIALS AND METHODS 

Cell culture and reagents 

Human cervical-carcinoma cells (HeLa) and Human em- 
bryonic kidney 293 (HEK293) cells were grown in 
Dulbecco's modified Eagle's medium (DMEM) supple- 
mented with 10% (v/v) fetal bovine serum and lOOU/ml 
penicillin and lOOmg/ml streptomycin and cultured at 
37°C and 5% C0 2 in a humidified incubator. 

For hypoxia treatments, cells were grown at 37°C in 
sealed chambers (Billups-Rothenberg) flushed with a 1% 
0 2 , 5% C0 2 , 94% N 2 gas mixture or in a Whitley 
hypoxystation (don Whitley Scientific, UK) set at 1% 
oxygen concentration. 

Dimethyloxalylglycine (DMOG, Frontier Scientific, 
CA, USA) was added to the indicated cultures at a 
500 uM final concentration. For the analysis of 
RUVBL2 and GYS1 expression, cDNA obtained from 
the following cell lines were also used: HepG2, HepaCl, 
NIH/3T3; N2a, HEK293 and A549. 

Analysis of gene profiling datasets 

The expression profiles corresponding to the indicated 
datasets and series were downloaded from Gene 
Expression Omnibus (GEO, http://www.ncbi.nlm.nih 
.gov/geo/) (28) database. In all the cases, untreated 
normoxic cells were used as reference. For each data set 
we calculated the mean of probe values in the biological 
replicates. All probes mapping to the locus of interest, 
except those with null values, were included in the 
analysis. Then, for each probe, the effect of hypoxia was 
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calculated as the logarithm of the ratio of the means of 
treated and control samples. Finally, individual log-ratio 
values were normalized by subtraction of the mean of all 
the log-ratios across the data set and division by their 
standard deviation. In the case multiple probes mapping 
to a given gene locus, the average of the log ratio was 
calculated. Information regarding GEO and probes ID 
can be found in the Supplementary Table SI. 

RNA extraction and quantitative PCR 

Total RNA was extracted and purified with the RNeasy 
Mini Kit (Qiagen). One microgram of RNA from each 
sample was reverse-transcribed to cDNA (Improm-II 
reverse transcriptase; Promega) and lul of cDNA 
samples were used as template for amplification reactions 
carried out with the LC Fast Start DNA master SYBR 
Green I kit (Roche Applied Science), following the manu- 
facturer's instructions. PCR amplifications were carried 
out in a Light Cycler System (Roche Applied Science), 
and data were analyzed with LightCycler software 3 
version 3.5.28 (Idaho Technology Inc.). For each 
sample, duplicate determinations were made, and the 
gene expression determined by the AACt method using 
(3-actin as reference gene. The primers used in this study 
are shown in the Supplementary Table S2. 

Plasmid construction 

Human genomic DNA extracted from HeLa cells was 
used as template for PCR amplification of GYS1 and 
RUVBL2 promoter regions. For reporter assays, PCR 
products were first cloned into pCR2.1-TOPO (invitrogen) 
and subsequently subcloned into the KpnI/XhoI restric- 
tion sites of pGL3-Basic. 

For the in vitro EBA, the putative insulators were 
cloned in the plasmid pELuc (29) and assayed their 
activity by transfection of the resulting constructs into 
HEK293 cells. This vector carries a CMV enhancer and 
a minimal promoter controlling the firefly luciferase report 
gene cassette with the polyadenylation site from SV40. All 
inserts were cloned into the Xhol site, between the 
enhancer and the CMV promoter (IN-position), and 
into Smal, upstream of the enhancer (OUT-position). 
The insertion in Smal site is a control for the assay that 
tests the potential silencer/repressor effects of our insert. 
For the in vivo EBA the putative insulator elements were 
cloned in the vector 48RCar (30) and injected resulting 
constructs into zebrafish embryos to assess their effect 
on transcription. In this plasmid the GFP gene is carried 
in the vector under the control of an actin promoter, that 
directs expression in heart and muscle. The insulator is 
cloned in the Kpnl site, between the actin promoter and 
the Enhancer 48, which targets the expression to the 
central nervous system (31). 

The identity of all constructs was verified by sequencing. 
All primer sequences are available in the Supplementary 
Table SI. 

Reporter assays 

Reporter assays were performed using the HeLa. Cells 
were seeded on six-well plates (3 x 10 5 cells/well) 6h 



prior to transfection. A 9 ug DNA mixture containing 
3 ug of the indicated reporter construct or empty 
plasmid and 0.6 ug of a plasmid encoding for renilla (sea 
pansy) luciferase under the control of a null promoter 
(Promega, Madison, WI, USA.) was used for transfection 
using the calcium phosphate method. On completing 16 h 
after transfection, cells where washed, replated in 24-well 
plates and incubated in normoxia, in the presence of 
DMOG or under hypoxia for additional 16 h. After treat- 
ments, the cells were lysed and the firefly and renilla 
luciferase activities of the lysate were determined using a 
dual-luciferase system (Promega, Madison, WI, USA). 
The firefly luciferase activity was normalized to that of 
renilla luciferase. 

In vitro EBA 

Cells were transiently transfected using Lipofectamine 
2000 (invitrogen) and OPTI-MEM® medium 
(Invitrogen) according to the supplier's instructions. 
Briefly, 1.8 xlO 5 HEK293 cells were seeded the day 
before transfection in 24-well plates. For each well, 
0.66 ug of the linearized reporter vector was transfected 
together with 0.14ug of pCMV-lacZ control plasmid (for 
normalization purposes). The pELuc-derived constructs 
were linearized prior to transfection, to avoid bidirectional 
enhancer activity, using the restiction enzyme Asp-718. 
This site is located 3' downstream from the 
polyadenylation signal. The pCMV-lacZ plasmid was 
also linearized using Seal. The cells were incubated 24 h 
with the transfection mixture and were thereafter lysed 
with Reporter Lysis Buffer (Promega). 

The luciferase activity was determinad using the 
Luciferase Assay Reagent (Promega) according to the 
manufacturer's specifications in a microplate luminometer 
(Orion, Berthold Detection Systems). The sample 
luciferase activity was corrected by the P-galactosidase 
activity in and the number of molecules (picomol) of the 
transfected plasmid construct (according to each plasmid 
size). Finally, activities were normalized as a fraction of 
the mean luciferase values obtained for the empty (pEluc) 
plasmid. 

Transgenesis and in vivo EBA 

Transgenesis and in vivo EBA were performed as described 
(30,31). For zebrafish transgenesis, the Tol2 transposon/ 
transposase method of transgenesis (32) was used with 
minor modifications. In total, 1 nl was injected in the 
cell of one-cell stage embryos containing 50ng/ul of 
transposase mRNA, 40ng/ul of phenol/chloroform 
purified DNA and 0.05% phenol red. 

Statistical analysis of data 

Statistical analysis of the experimental data was per- 
formed with the R software package [(33), http://www 
.R-project.org/]. The statistical tests applied to each data 
set are indicated in the figure legends. We adopted the 
following code to indicate the magnitude of P-values 
throughout the manuscript figures as: ***p = [0, 0.001]; 
**P = [0.001, 0.01] and*i> = [0.01, 0.05]. 
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RESULTS 

Differential regulation of GYS1 and RUVBL2 by hypoxia 

The muscle glycogen synthase gene, GYS1, is regulated by 
HIF as part of the cellular metabolic reprogramming 
required for the adaptation to hypoxia (27). The regula- 
tion of human GYS1 by hypoxia is mediated by a func- 
tional RCGTG element located 255 bp upstream its 
transcription start site (TSS) (27). Very close to GYS1, 
but encoded by the opposite DNA strand, is located the 
RUVBL2 gene (Figure 1A). The relative position, 
intergenic distance and orientation of these two genes 
are conserved across mammals (data not shown). The 
TSS of RUVBL2 is located at 288 bp of the HRE 
driving GYS1 expression in response to hypoxia 
(location of the HRE is shown by a black box in the 
'blaf track, Figure 1A), raising the possibility of a 
coordinated regulation of these two genes by HIF. In 
fact, the intergenic region between GYS1 and RUVBL2 
can be considered a bidirectional promoter (Figure 1A, 
'Elnitski bidirectional promoters' prediction track). To 
study this possibility, we analyzed publicly available 
gene expression profiles of cells exposed to hypoxia and 
found that, whereas GYS1 mRNA levels were induced by 
hypoxia in most of the profiles, the expression of RUVBL2 
remained constant or was even repressed under low 
oxygen tension (Figure IB). To confirm these results we 
exposed myotubes to hypoxia and determined its effect on 
GYS1 and RUVBL2 mRNA levels. As shown in 
Figure 1C, GYS1 mRNA level significantly increased in 
response to hypoxia, in agreement with published results 
(27), but the treatment did not induce RUVBL2 mRNA. 
In order to rule out a cell-type specific effect, we 
determined the effect of hypoxia on the expression of 
GYS1 and RUVBL2 in a variety of cell types 
(Figure ID). These analyses confirmed that, whereas 
GYS1 expression was increased by hypoxia in virtually 
all cell lines studied, RUVBL2 remained largely unaffect- 
ed. In agreement with this conclusion, a meta-analysis of 
gene expression profile experiments suggested that, unlike 
GYS1, RUVBL2 is not significantly modulated by hypoxia 
(10). These results indicated that hypoxia specifically 
affects GYS1, but not RUVBL2 transcription. To 
confirm this possibility, we investigated the effect of 
hypoxia on a set of reporter constructs derived from this 
locus. As shown in Figure IE, a reporter construct con- 
taining the region flanking GYS1 gene and including the 
HRE (region cG spanning residues +84 to —429 relative to 
GYS1 TSS, Figure 1A), was robustly induced by hypoxia 
and the hypoxia mimetic DMOG in HeLa cells, as 
expected from the presence of the evolutionarily conserved 
block containing the HRE within this cG region. In 
contrast, the analogous region upstream of RUVBL2 
(cR, spanning -396 to +12 relative to RUVBL2 TSS) 
did not respond to HIF activation in spite of including 
the same HRE-containing region (Figure IE). 
Importantly, this same result was obtained when the 
whole intergenic region, maintaining the original 
genomic structure, was used to drive luciferase expression 
from the GYS1 (cRcG) or the RUVBL2 (cGcR) pro- 
moters (Figure IE). 



Collectively, these results indicate that the HRE located 
between GYS1 and RUVBL2 genes selectively drives the 
transcription of the former in response to hypoxia. 

The lack of RUVBL2 induction by hypoxia is not due 
to gene silencing 

The accessibility of promoter regions is one that the mech- 
anisms that restricts the activity of a TF toward specific 
genes. This can be achieved by methylation of promoter 
regions and/or by altering chromatin compactness 
through histone modification. In fact, it has been 
recently found that basal promoter activity determines 
cell type-specific HIF transcription (8). Thus, we 
investigated whether a lack of RUVBL2 transcriptional 
activity could explain the observed selectivity of hypoxia 
within this locus. As shown in Figure 2A, the normoxic 
levels of RUVBL2 and GYS1 mRNAs were of similar 
magnitude, at least for the set of cell lines included in 
our study. Moreover, the differential effect of hypoxia 
(Figure ID) was observed even for cell lines, such as 
HeLa or A549, in which the relative basal level of 
RUVBL2 mRNA was much higher than that of GYS1. 
On the other hand, the genomic region adjacent to the 
RUVBL2 gene (cR) showed a strong basal promoter 
activity when assayed in HeLa cells (Figure 2B). In fact, 
the promoter activity of the cR region was significantly 
higher than that of the cG region (Figure 2B). In agree- 
ment with our results, published data of genome-wide 
RNA polymerase II binding shows a similar signal in 
the GYS1 and RUVBL2 promoter regions in a wide 
range of cell types (Supplementary Figure SI). 

Altogether, these results indicate that RUVBL2 and 
GYS1 are transcribed under normoxia to a similar 
extent and thus, we discarded RUVBL2 promoter accessi- 
bility as a potential explanation for the differential regu- 
lation observed under hypoxia. 

An inhibitory region prevents RUVBL2 induction 
by the HRE 

To gain insight into the molecular mechanism responsible 
for the differential regulation of RUVBL2 and GYS1 by 
hypoxia, we generated a set of reporter constructs contain- 
ing different deletions and rearrangements of the genomic 
region between GYS1 and RUVBL2 genes (Figure 3A and 
diagrams on the left of Figure 3B). Based on the evolu- 
tionary conservation (Phas Cons elements, see Figure 1A), 
we differentiated five blocks within this region (Figures 1A 
and 3): proximal GYS1 (pG), upstream GYS1 (uG), 
HRE-containing block (HRE), upstream RUVBL2 (uR) 
and proximal RUVBL2 (pR). As shown in Figure 3, 
analysis of the transcriptional activity of these reporter 
constructs in HeLa cells showed that neither pG nor pR 
proximal regions responded to HIF activation. However, 
combination of these proximal regions with the HRE 
block, regardless of the orientation of the latter, resulted 
in constructs (HREF_pG, HRER_pG, HREF_pR and 
HRER_pR) that were robustly induced by the hypoxia 
mimetic DMOG. The induction of these constructs was 
of similar magnitude to that observed for the complete 
GYS1 construct (cG), suggesting that upstream regions 
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(uG and uR) are not required for the activity of the 
HRE-containing block. These data ruled out a collabor- 
ation between HIF and other transcription factors binding 
to pG or uG as a potential explanation for the differential 
induction of GYS1 and RUVBL2 by hypoxia. On the 
other hand, comparison of response of the cR and 
HREF_pR constructs, hinted the existence of an inhibi- 
tory region within the uR region as removal of this region 
significantly increased the response to HIF activation. In 
agreement, the up-regulation of the HREF_pG construct 
was blunted by the inclusion of the upstream region of 
RUVBL2 between the HRE and proximal GYS1 block 
(Figure 3, HREF_pG and HREF uR pG). 

In summary, these results suggest that the genomic 
region from -202 to -30, relative to the RUVBL2 TSS, 
(uR) prevents the induction of RUVBL2 by HIF. 

The GYS1/RUVBL2 intergenic region contains an 
enhancer blocking element 

Insulators are DNA elements that can prevent the promis- 
cuous effect of enhancers or silencers, restricting their 
interactions with promoters. Our results suggested that 
the uR region prevents the interaction of the HRE 
over the RUVBL2 promoter, a function that is compatible 
with the enhancer blocking type activity commonly 
associated to insulators. To test this possibility, we 
studied the ability of different sequences, derived from 
the RUVBL2/GYS1 intergenic region, to interfere with 
the activity of a heterologous enhancer/promoter pair in 
a standard EBA, in HEK293 cells (29). Each of the 
genomic fragments under study was cloned either, 
between the enhancer and promoter (IN position) or 



upstream and the enhancer (OUT position) (Figure 4A). 
As a positive control, we also included the boxes 11 /III (E 
II/III) derived from the classical chicken 5'HS4 P-globin 
insulator element (34), known to bind the nuclear factor 
CTCF and responsible for the enhancer blocking effect of 
the 5'HS4 element (35). The enhancer blocking activity of 
these sequences, represented as fold reduction of the 
activity of a vector lacking insert, is represented in 
Figure 4B. In agreement with previous reports (35), the 
E II/III boxes, but not a mutant form, interfered with the 
activity of the CM V enhancer when inserted into the Xhol 
site (Figure 4B, white bars). Enhancer blockers work only 
in the 'IN' configuration as demonstrated by the lack of 
activity of the EII/III constructs in the 'OUT' configur- 
ation (Figure 4B, black bars). Significantly, all the con- 
structs containing the uR in the 'IN' configuration, except 
1R, showed a significant reduction in the enhancer activity 
(Figure 4B, white bars corresponding to constructs IF, 
2F, 2R, 3F and 3R). Importantly, these same regions 
had no significant inhibitory effect when cloned in the 
'OUT' configuration (Figure 4B, black bars corresponding 
to constructs IF, 2F, 2R, 3F and 3R). In contrast to these 
results, neither the HRE nor uG regions were able to 
suppress the activity of the CMV enhancer (Figure 4B, 
constructs 4F, 4R, 5F, 5R, 6F, 6R and Supplementary 
Figure S2). 

The EBA results strongly support the existence of an 
enhancer blocking element within the genomic region 
from -202 to -30 relative to the RUVBL2 TSS 
(Figure 4B, constructs 3F and 3R). Thus, we decided to 
functionally validate the presence of an insulator in this 
region in vivo by means of an independent assay using 



Figure 1. Continued 

df = 7, P = 0.002964) from the value of one (no induction). (E) HeLa cells were transiently transfected with reporter constructs containing the 
indicated region (cG, cR, cGcR or cRcG, see Figure 1A) upstream a firefly luciferase gene and exposed to normoxia, hypoxia or the hypoxia mimetic 
DMOG for 12 h. The graphs represent the corrected luciferase activity values of each construct in cells exposed to hypoxia/DMOG and represented 
as fold change over the activity obtained in normoxic cells. Bars represent the average of values obtained in three independent experiments and error 
bars their standard deviation. Statistically significant differences between pairs of constructs are indicated by asterisks (one-way ANOVA, 
F9,2o = 36.704, P = 1.6 x 10~'°, followed by Tukey's multiple comparison test). 
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Figure 3. Region uR prevents the increase in transcription mediated by the HRE. (A) Schematic diagram of the GYS1/RUVBL2 intergenic showing 
the different elements included in the reporter constructs. (B) HeLa cells were transfected with the indicated constructs and treated with the hypoxia 
mimetic DMOG or left untreated for 12 h. The graph represents the corrected luciferase activity values of each construct in treated cells represented 
as fold change over the activity obtained in control (normoxic) cells. Bars represent the average of values obtained in three independent experiments 
and error bars, their standard deviation. The diagrams on the left are a schematic representation of the genomic elements included in the reporter 
construct and their orientation relative to the luciferase gene. Statistically significant differences between pairs of constructs within the groups 
containing pG or pR are indicated by asterisks (one-way ANOVA followed by Tukey's multiple comparison test. pG-containing constructs 
F 4i i 0 = 13.608, P = 0.0004705; pR-containing constructs F 4JQ = 7.518, P = 0.004598). 



transgenic zebrafish. To this end, we used a reporter con- 
struct in which EGFP expression is under the control of 
the cardiac actin promoter from Xenopus laevis and the 
Z48 enhancer, which drives transgene expression in the 
midbrain (Figure 5 A) (31). Transgenic zebrafish embryos 
injected with this construct showed EGFP expression in 
the heart and in the developing somites (due to the cardiac 
actin promoter), as well as in the CNS (Figure 5B and C, 
control), due to the Z48 enhancer, as reported earlier 
(30,31). Micro-injection of a construct containing the uR 
region cloned between the enhancer and the promoter, 
regardless of its orientation, resulted in a strong attenu- 
ation of the CNS signal, whereas retaining the EGFP 
signal in the developing somites (Figure 5C, uR_F and 
uR_R). Analysis of the somites/CNS EGFP-mediated 
fluorescence signal ratio in 62 independent transgenic 
zebrafish lines demonstrated an enhancer blocking 
activity associated with the uR sequence (Figure 5B, 
uR_F and uR_R). Importantly, in agreement with the 
in vitro EBA assays (Figure 4), the effect was specific for 
the uR region as cloning of the uG sequence between the 
Z48 enhancer and the actin promoter did not interfere 
with the EGFP expression in the CNS (Figure 5B and 
C, uG_F and uG_R). 



Altogether, these results demonstrate the existence of a 
powerful enhancer blocking element in the uR region that 
could explain the differential regulation of GYSlj 
RUVBL2 by HIF. The activity of this insulator element 
is observed when assayed in a distant heterologous system 
(zebrafish) suggesting that its function and potential trans- 
acting factors have been evolutionary conserved in 
vertebrates. 



DISCUSSION 

The recognition of short DNA motifs by transcription 
factors is a key step in the regulation of transcription. 
However, the low information content of most of the 
TFBS predicts a promiscuous binding that is in contrast 
with the observed specificity. This apparent paradox raises 
the question of how transcription factor selectivity is 
achieved. Although several mechanisms including chro- 
matin accessibility and TF cooperation, have been 
shown to restrict the target space for a given TF, a 
complete explanation is still lacking in most of the cases. 
In this scenario, HIF transcription factors are not an ex- 
ception. HIF heterodimers bind to the RCGTG motif 
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Figure 4. The uR region contains an enhancer blocking element. (A) Schematic diagrams of the different elements within the GYS1/RUVBL2 
intergenic region (left diagram) and the EBA vector showing the 'IN' (middle diagram) and 'OUT' (right diagram) cloning sites. (B) Each of the 
indicated elements (diagrams on the left under 'INSERT') derived from the GYS1/RUVBL2 intergenic region, or from the (3-globin gene (E II/III and 
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vector to generate the indicated constructs ('CONST.'). The HEK293 cells were transiently transfected with each one of these constructs, an empty 
EBA vector lacking insert (pEluc) or an empty EBA vector lacking the enhancer (pAEluc). On completing 24 h after transfection, the cells were 
processed to determine the transcriptional activity of the constructs. The figure represents the corrected luciferase activity in each sample and is 
expressed as fold reduction of the activity observed in cells transfected with the empty vector, pEluc (29). Bars represent the average of values 
obtained in three independent experiments and errors bars, their standard deviation. Asterisks indicate constructs whose activity was significantly 
different from that observed for the control, pEluc (one-way ANOVA, ^29,237 = 131.10, P = 2.2 x 1CT 16 , followed by Tukey's multiple comparison 
test). 



(4,6,9), which is present in almost every gene in the human 
genome, yet hypoxia results in the regulation of a few 
hundred genes only and, accordingly, HIF only binds to 
a subset of the potential binding sites (5,6,9,10). It has 
been recently demonstrated that HIF binds preferentially 
to RCGTG motifs present in the promoters of genes 
actively transcribed under normoxic conditions (8). 
Although this restriction results in a large reduction in 
the number of potential targets, it does not fully explain 
HIF selectivity as many genes transcribed under normoxia 
and containing RCGTG motifs are not induced by 
hypoxia. With the aim of gaining insight into the mech- 
anisms of HIF target selectivity, we investigated the 
GYS1 /RUVBL2 genes as a particularly striking example 
of differential regulation. Given their close proximity and 
the location of the functional HRE between both genes, 
we had expected their coordinated regulation by HIF. 
However, our results showed that while GYS1 was 



induced by hypoxia, RUVBL2 levels remained unchanged 
(Figure 1). The basal level of mRNA, RNA polymerase II 
binding and promoter activity (Figure 2 and 
Supplementary Figure SI) ruled out the accessibility of 
chromatin as an explanation for the lack of induction of 
RUVBL2 by hypoxia. Detailed analysis of the promoter 
activity of different fragments derived from the GYSlj 
RUVBL2 intergenic region revealed the existence of an 
inhibitory region between the HRE and the RUVBL2 
gene (Figure 3). Finally, specific EBAs performed 
in vitro (Figure 4), as well as in vivo (Figure 5) 
demonstrated the existence of an insulator element that 
could explain the lack of effect of the HRE over the 
RUVBL2 promoter and thus, the differential regulation 
of GYS1 and RUVBL2 genes by hypoxia. 

Insulators have been shown to play a key role in the 
differential patterns of gene expression during develop- 
ment and cell-lineage specification (23). However, to our 
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Figure 5. The insulator element derived from RUVBL2 genomic region blocks enhancer activity in vivo. (A) Schematic representation of the ZED 
construct. (B and C) zebrafish embryos were micro-injected with an empty ZED vector (control), lacking insert between the enhancer and promoter, 
or constructs containing the RUVBL2-derived enhancer blocking element (uR) or the corresponding region upstream the GYS1 gene (uG). Both 
elements were cloned in the forward (uR_F, uG_F) or reverse (uR_R, uG_R) orientation. After 36 h after micro-injection, the GFP signal in CNS 
and somites was determined by fluorescence microscopy and quantified with LaserPix (Bio-Rad) image analysis software, as reported before (30). 
(B) A representative image of each group of animals is shown. Arrowheads indicate the location of the CNS. (C) The graph represents the ratio 
between the GFP signals in somites and CNS for each of the transgenic fish. Horizontal line represents the median of the ratios for each group of 
animals. Asterisks indicate set of values that were significantly different from those obtained for the empty construct Kruskal-Wallis, 
chi-squared = 89.2081, df = 5, P<2.2e-16; P-values for comparisons were calculated by Wilcoxon rank and corrected for multiple comparison 
using Bonferoni's method. 



knowledge, this is the first example of an enhancer 
blocking element contributing to the selectivity of a 
transcription factor acutely induced by environmental 
factors, such as HIF. Although further work is required 
to determine whether this is a general mechanism 
contributing to HIF specificity, our working hypothesis 
is that this mechanism could be particularly relevant in 
cases of bidirectional promoters, such as the one described 
herein, where chromatin accessibility and DNA methyla- 
tion is likely to be similar across the intergenic region. In 
support to this possibility, we have found at least one 
further example of bidirectional promoters differentially 
regulated by hypoxia, BCKDHA and EXOSC5. 
Preliminary results indicate that, in spite of similar basal 
transcription, BCKDHA, but not EXOSC5, is induced by 
hypoxia (Supplementary Figure S3). 

Another issue raised by our results relates to the mo- 
lecular identity of the insulator element located upstream 
of the RUVBL2 gene. In vertebrates, several regulatory 
elements including CTCF binding motifs (36-38), 



repetitive elements, [such as ALUs (39), SINE B2 (29) 
and SINE Bl (30)] and scaffold/matrix-attachment 
regions [S/MARs; (40,41)], have been shown to function 
as insulators (25,42). Among them, the most widespread 
and well characterized are CTCF-binding elements (43). 
We have not found obvious CTCF binding motif within 
the —202 to —30 region upstream of RUVBL2, and pub- 
lished ChIP data shows CTCF binding to the GYS1 TSS 
region, but not to the region between the HRE and 
RUVBL2 gene (Supplementary Figure S4A, CTCF ChIP 
signal track). On the other hand, there are no repetitive 
elements within this genomic region (Supplementary 
Figure S4A, repeats tracks) arguing against the existence 
of a SINE element that could explain the observed 
enhancer blocking activity. Finally, although we found 
no locally high proportion of A/T nucleotides, typically 
associated with S/MARs elements (41) in the region 
upstream of RUVBL2 (Supplementary Figure S4A, GC 
percent track), computer prediction of S/MAR sites, 
found a significant signal in the upstream region of 
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RUVBL2 (Supplementary Figure S4B). Thus, the insula- 
tor activity described herein could be mediated by S/MAR 
elements. However, as this evidence has been obtained 
through in silico approaches, further work is required to 
identify the minimal region required for the enhancer 
blocking function, including its molecular identity, 
underlying mechanism and the associated potential 
trans-acting factors. This is yet another perfect example 
demonstrating the diversity of mechanisms, most of them 
not known to date, that cells are using to organize func- 
tional insulator elements (25). 

The mechanism by which enhancer blockers prevent the 
activity of upstream enhancers is unclear, although our 
current understanding is that insulators will probably 
not be using unique mechanisms but, rather, adaptations 
of pre-existing ones already in place for the normal regu- 
lation of gene expression (23). Several of the proposed 
models invoke the generation of chromatin loops that seg- 
regate enhancer-sensitive and resistant promoters in 
distinct domains (44,45). In the case described herein, 
the model is further complicated by the short distance 
existing between the ra-elements involved. A piece of in- 
formation that could shed light into the mechanism is the 
intriguing observation that the inclusion of the uG region 
seems to abolish the enhancer blocking effect of uR in 
EBA (Figure 4B, compare constructs 1R/1F with 3R/ 
3F). Although we cannot currently explain this behavior, 
it would suggest that the EBA of uR can be modulated by 
elements located in its proximity. In this regard, it has 
been previously found that the EBA of the 'gypsy' 
element is affected by the number of copies of this 
element. When two copies, instead of one, are located 
between the enhancer and promoter, its blocking effect 
is abolished (46). However, the enhancer blocking effect 
is restored by the insertion of a third copy (47), depending 
on the order and distance of the insulator elements, 
indicating that complex protein-protein interactions are 
responsible for these unexpected effects and underlying a 
major role of insulators in whole nuclear genome organ- 
ization (24,48). Thus, it is possible that the impairment of 
the EBA of uR by uG can be reverted by other cis- 
elements, present in the native genomic context, but not 
included in this set of constructs. A further possibility is 
that the HRE-containing block could enhance the tran- 
scription from the minimal CMV promoter. If this is the 
case, the enhancer blocking element within uR would be 
located upstream of the HRE element in the construct 1R 
and would be thus, unable to prevent its action upon the 
minimal CMV promoter. Regardless of the specific mech- 
anism by which uR exerts its effect, our data clearly dem- 
onstrates that it contains an EBA. 

Reprogramming of cellular metabolism, in particular 
glucose metabolism, is central in the cellular adaptation 
to hypoxia. The hypoxic induction of GYS1, encoding for 
an isoform of glycogen synthase, is part of this 
reprogramming (27). On the other hand, the existence of 
an enhancer blocking element between the HRE and the 
RUVBL2 promoter raises the question of why RUVBL2 
expression has to be shielded from the HIF-mediated in- 
duction. RUVBL2 gene encodes for Reptin, an AAA+ 
ATPase that forms part of chromatin remodeling 



complexes. Interestingly, it has been recently shown that 
hypoxia leads to reptin methylation and that, upon this 
post-translational modification, it is able to repress 
HIF-mediated transcription (49). Thus, it is plausible 
that RUVBL2 expression might not be induced by 
hypoxia so as not to upset the balance between 
methylated/unmethylated reptin that could lead to prema- 
ture termination of HIF-mediated transcription. However, 
this hypothesis does not provide an explanation for the 
close proximity of these genes. The conservation of the 
GYS1/RUVBL2 genomic arrangement from opossums to 
humans suggests a selective pressure to maintain both 
genes in close proximity. Previous studies have shown 
that a substantial proportion of mammalian genes is 
arranged in a divergent head-to-head structure and 
controlled by bidirectional promoters (50), so that the 
pair of genes tend to be co-expressed (50,51). The need 
for co-regulation of the pair of genes under the control 
of a bidirectional promoter could explain the selective 
pressure that keeps them in close proximity. However, in 
the case of GYS1 / RUVBL2, co-regulation of both genes is 
an unlikely reason for their close proximity as the exist- 
ence of the insulator element would prevent the action of 
flanking c/.v-elements on the opposite promoter. 

In summary, we have identified an insulator, acting as a 
functional enhancer blocking element, that explains the 
differential response of GYS1 and RUVBL2 genes to 
hypoxia. To our knowledge, this is the first report 
describing a role for this type of genetic elements in 
dictating the specificity of acutely induced transcription 
factors in response to environmental, as opposed to devel- 
opmental, signals. Importantly, the generalization of this 
model adds to the arsenal of strategies, including chroma- 
tin accessibility and combinatorial assembly of TFs, that 
are employed by eukaryotes to ensure a highly specific 
gene expression based on an otherwise promiscuous set 
of c«-regulatory elements. 

SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online: 
Supplementary Tables 1 and 2 and Supplementary 
Figures 1^1, Supplementary Reference [54]. 

ACKNOWLEDGEMENTS 

The data used for Supplementary Figures SI and S4 were 
generated and analyzed by the labs of Michael Snyder, 
Mark Gerstein and Sherman Weissman at Yale 
University; Peggy Farnham at UC Davis and Kevin 
Struhl at Harvard (http://genome.ucsc.edu/cgi-bin/ 
hgTrackUi?hgsid = 19329803 l&c = chrl&g = wgEncode 
YaleChlPseq) as members of the the ENCODE Project 
consortium (52). Luis del Peso and Lluis Montoliu 
contributed equally to this work. 

FUNDING 

The Ministerio de Ciencia e Innovation (Spanish Ministry 
of Science and Innovation, MICINN) grant numbers 



7926 Nucleic Acids Research, 2012, Vol. 40, No. 5 



SAF2008-03147 (to L.del P.), BIO2009-1297 (to L.M.), 
BFU2010-14839 (to J.L.G.-S.), CSD2007-00008 (to 
J.L.G.-S.); Comunidad Autonoma de Madrid grant 
number S-SAL-0311_2006 (to L.del P.); the 7th 
Research Framework Programme of the European 
Union grant number METOXIA project ref. 
HEALTH-F2-2009-222741 (to L.del P.); Junta de 
Andalucia grant number CVI-3488 (to J.L.G.-S.); 
CIBERER (ISCIII) (to E.M.). Funding for open access 
charge: SAF2008-03147 (to L.del P.). 

Conflict of interest statement. None declared. 



REFERENCES 

1. Wang,G.L., Jiang.B.H., Rue,E.A. and Semenza,G.L. (1995) 
Hypoxia-inducible factor 1 is a basic-helix-loop-helix-PAS 
heterodimer regulated by cellular 02 tension. 

Proc. Natl Acad. Sci. USA, 92, 5510-5514. 

2. Salceda.S. and CaroJ. (1997) Hypoxia-inducible factor lalpha 
(HIF-lalpha) protein is rapidly degraded by the 
ubiquitin-proteasome system under normoxic conditions. Its 
stabilization by hypoxia depends on redox-induced changes. 

J. Biol. Chem., 272, 22642-22647. 

3. Jiang,B.H., ZhengJ.Z., Leung,S.W., Roe,R. and Semenza,G.L. 
(1997) Transactivation and inhibitory domains of 
hypoxia-inducible factor lalpha. Modulation of transcriptional 
activity by oxygen tension. /. Biol. Chem., 272, 19253-19260. 

4. Wang,G.L. and Semenza.G.L. (1993) Characterization of 
hypoxia-inducible factor 1 and regulation of DNA binding 
activity by hypoxia. J. Biol. Chem., 268, 21513-21518. 

5. Mole,D.R., Blancher.C, Copley.R.R., Pollard,P.J., GleadleJ.M., 
Ragoussis,J. and Ratcliffe,P.J. (2009) Genome-wide association of 
hypoxia-inducible factor (HIF)- lalpha and HIF-2alpha DNA 
binding with expression profiling of hypoxia-inducible transcripts. 
J. Biol. Chem., 284, 16767-16775. 

6. Xia,X., Lemieux,M.E., Li,W., CarrollJ.S., Brown.M., Liu,X.S. 
and Kung,A.L. (2009) Integrative analysis of HIF binding and 
transactivation reveals its role in maintaining histone methylation 
homeostasis. Proc. Natl Acad. Sci. USA, 106, 4260-4265. 

7. Wenger.R.H., Stiehl,D.P. and Camenisch,G. (2005) Integration of 
oxygen signaling at the consensus HRE. Sci. STKE, 2005, rel2. 

8. Xia,X. and Kung,A.L. (2009) Preferential binding of HIF-1 to 
transcriptionally active loci determines cell-type specific response 
to hypoxia. Genome Biol., 10, R113. 

9. SchodelJ., Oikonomopoulos.S., RagoussisJ., Pugh,C.W., 
Ratcliffe,P.J. and Mole,D.R. (2011) High-resolution genome-wide 
mapping of HIF-binding sites by ChlP-seq. Blood, 117, 
e207-e217. 

10. Ortiz-Barahona,A., Villar,D., Pescador,N., AmigoJ. and del 
Peso,L. (2010) Genome-wide identification of hypoxia-inducible 
factor binding sites and target genes by a probabilistic model 
integrating transcription-profiling data and in silico binding site 
prediction. Nucleic Acids Res., 38, 2332-2345. 

11. Benita,Y., Kikuchi,H., Smith,A.D., Zhang,M.Q., Chung,D.C. and 
Xavier,R.J. (2009) An integrative genomics approach identifies 
Hypoxia Inducible Factor-1 (HIF-l)-target genes that form the 
core response to hypoxia. Nucleic Acids Res., 37, 4587-4602. 

12. Pan,Y., Tsai,C.J., Ma,B. and Nussinov,R. (2010) Mechanisms of 
transcription factor selectivity. Trends Genet., 26, 75-83. 

13. Ogata,K., Sato,K. and Tahirov,T.H. (2003) Eukaryotic 
transcriptional regulatory complexes: cooperativity from near and 
afar. Curr. Opin. Struct. Biol, 13, 40—48. 

14. Panne, D., Maniatis,T. and Harrison,S.C. (2007) An atomic model 
of the interferon-beta enhanceosome. Cell, 129, 1111-1123. 

15. Zhang.W., Tsuchiya,T. and Yasukochi,Y. (1999) Transitional 
change in interaction between HIF-1 and HNF-4 in response to 
hypoxia. J. Hum. Genet., 44, 293-299. 

16. Yamashita,K., Discher,D.J., Hu,J., Bishopric,N.H. and 
Webster.K.A. (2001) Molecular regulation of the endothelin-1 
gene by hypoxia. Contributions of hypoxia-inducible factor-1, 



activator protein- 1, GATA-2, AND p300/CBP. J. Biol. Chem., 
276, 12645-12653. 

17. Aprelikova,0., Wood,M., Tackett,S., Chandramouli.GV. and 
BarrettJ.C. (2006) Role of ETS transcription factors in the 
hypoxia-inducible factor-2 target gene selection. Cancer Res., 66, 
5641-5647. 

18. Hu,C.J., Sataur,A., Wang.L., Chen,H. and Simon.M.C. (2007) 
The N-terminal transactivation domain confers target gene 
specificity of hypoxia-inducible factors HIF-lalpha and 
HIF-2alpha. Mol. Biol. Cell, 18, 4528^1542. 

19. Lupien,M., EeckhouteJ., Meyer,C.A., Wang,Q., Zhang,Y., Li,W., 
CarrollJ.S., Liu,X.S. and Brown,M. (2008) FoxAl translates 
epigenetic signatures into enhancer-driven lineage-specific 
transcription. Cell, 132, 958-970. 

20. Iguchi-Ariga.S.M. and Schaffner,W. (1989) CpG methylation of 
the cAMP-responsive enhancer/promoter sequence TGACGTCA 
abolishes specific factor binding as well as transcriptional 
activation. Genes Dev., 3, 612-619. 

21. Perini,G, Diolaiti.D., Porro,A. and Delia Valle.G. (2005) In vivo 
transcriptional regulation of N-Myc target genes is controlled by 
E-box methylation. Proc. Natl Acad. Sci. USA, 102, 12117-12122. 

22. Wenger.R.H., KvietikovaJ., Rolfs,A., Camenisch,G. and 
Gassmann,M. (1998) Oxygen-regulated erythropoietin gene 
expression is dependent on a CpG methylation-free 
hypoxia-inducible factor-1 DNA-binding site. Eur. J. Biochem., 
253, 771-777. 

23. Gaszner,M. and Felsenfeld,G (2006) Insulators: exploiting 
transcriptional and epigenetic mechanisms. Nat. Rev. Genet., 7, 
703-713. 

24. Capelson,M. and Corces,V.G. (2004) Boundary elements and 
nuclear organization. Biol. Cell, 96, 617-629. 

25. Molto.E., Fernandez,A. and Montoliu,L. (2009) Boundaries in 
vertebrate genomes: different solutions to adequately insulate gene 
expression domains. Brief. Fund. Genomic Proteomic, 8, 283-296. 

26. Burgess-Beusse,B., Farreil.C, Gaszner,M., Litt,M., Mutskov,V., 
Recillas-Targa,F., Simpson,M., West,A. and Felsenfeld,G. (2002) 
The insulation of genes from external enhancers and silencing 
chromatin. Proc. Natl Acad. Sci. USA, 99, 16433-16437. 

27. Pescador,N., Villar,D., Cifuentes,D., Garcia-Rocha,M., 
Ortiz-Barahona,A., Vazquez,S., Ordonez,A., Cuevas,Y., 
Saez-Morales,D., Garcia-Bermejo,M.L. et al. (2010) Hypoxia 
promotes glycogen accumulation through hypoxia inducible factor 
(HIF)-mediated induction of glycogen synthase 1. PLoS One, 5, 
e9644. 

28. Barrett,T., Troup.D.B., Wilhite,S.E., Ledoux.P., Rudnev.D., 
Evangelista,C, KimJ.F., Soboleva,A., Tomashevsky,M., 
Marshall,K.A. et al. (2009) NCBI GEO: archive for 
high-throughput functional genomic data. Nucleic Acids Res., 37, 
D885-D890. 

29. Lunyak,V.V., Prefontaine,G.G., Nunez,E., Cramer,T., Ju,B.G., 
Ohgi,K.A., Hutt,K., Roy,R., Garcia-Diaz,A., Zhu,X. et al. (2007) 
Developmentally regulated activation of a SINE B2 repeat as a 
domain boundary in organogenesis. Science, 317, 248-251. 

30. Roman.A.C, Gonzalez-Rico, F J., Molto,E., Hernando,H., 
Neto.A., Vicente-Garcia,C, Ballestar.E., Gomez-SkarmetaJ.L., 
Vavrova-Anderson.J., White,R.J. et al. (2011) Dioxin receptor 
and SLUG transcription factors regulate the insulator activity of 
Bl SINE retrotransposons via an RNA polymerase switch. 
Genome Res., 21, 422^132. 

31. BessaJ., TenaJ.J., de la Calle-Mustienes,E., Fernandez-Minan.A., 
Naranjo,S., Fernandez.A., Montoliu,L., Akalin,A., Lenhard,B., 
Casares,F. et al. (2009) Zebrafish enhancer detection (ZED) 
vector: a new tool to facilitate transgenesis and the functional 
analysis of cis-regulatory regions in zebrafish. Dev. Dyn., 238. 
2409-2417. 

32. Kawakami.K. (2004) Transgenesis and gene trap methods in 
zebrafish by using the Tol2 transposable element. Methods Cell 
Biol, 11, 201-222. 

33. R Development Core Team. (2008) R: A Language and 
Environment for Statistical Computing. Vienna, Austria. 

34. Recillas-Targa,F., Bell,A.C. and Felsenfeld.G. (1999) Positional 
enhancer-blocking activity of the chicken beta-globin insulator in 
transiently transfected cells. Proc. Natl Acad. Sci. USA, 96, 
14354-14359. 



Nucleic Acids Research, 2012, Vol. 40, No. 5 1927 



35. Recillas-Targa,F., Pikaart.M.J., Burgess-Beusse.B., Bell,A.C, 
Litt,M.D., West,A.G., Gaszner,M. and Felsenfeld.G. (2002) 
Position-effect protection and enhancer blocking by the 
chicken beta-globin insulator are separable activities. 

Proc. Natl Acad. Sci. USA, 99, 6883-6888. 

36. Beh\A.C, West,A.G. and Felsenfeld.G. (1999) The protein 
CTCF is required for the enhancer blocking activity of vertebrate 
insulators. Cell, 98, 387-396. 

37. Furlan-Magaril,M., Rebollar,E., Guerrero,G, Fernandez,A., 
Molto.E., Gonzalez-Buendia,E., Cantero,M., Montoliu,L. and 
Recillas-Targa,F. (2011) An insulator embedded in the chicken 
alpha-globin locus regulates chromatin domain configuration and 
differential gene expression. Nucleic Acids Res., 39, 89-103. 

38. Giraldo,P., Martinez,A., Regales,L., Lavado,A., Garcia-Diaz,A., 
Alonso,A., Busturia,A. and Montoliu.L. (2003) Functional 
dissection of the mouse tyrosinase locus control region identifies 
a new putative boundary activity. Nucleic Acids Res., 31, 
6290-6305. 

39. Willoughby.D.A., Vilalta,A. and Oshima,R.G. (2000) An Alu 
element from the K18 gene confers position-independent 
expression in transgenic mice. /. Biol. Chem., 275, 759-768. 

40. McKnight,R.A., Shamay,A., Sankaran,L., Wall,R.J. and 
Hennighausen,L. (1992) Matrix-attachment regions can impart 
position-independent regulation of a tissue-specific gene in 
transgenic mice. Proc. Natl Acad. Sci. USA, 89, 6943-6947. 

41. Millot,B., Montoliu,L., Fontaine,M.L., Mata,T. and Devinoy,E. 
(2003) Hormone-induced modifications of the chromatin structure 
surrounding upstream regulatory regions conserved between the 
mouse and rabbit whey acidic protein genes. Biochem. J., 372, 
41-52. 

42. Lunyak.V.V. (2008) Boundaries. Boundaries . . . Boundaries???. 
Curr. Opin. Cell Biol, 20, 281-287. 

43. West.A.G., Gaszner,M. and Felsenfeld.G. (2002) Insulators: many 
functions, many mechanisms. Genes Dev., 16, 271-288. 

44. Farrell,C.M., West,A.G. and Felsenfeld,G. (2002) Conserved 
CTCF insulator elements flank the mouse and human beta-globin 
loci. Mol. Cell Biol, 22, 3820-3831. 

45. Byrd.K. and Corces.V.G. (2003) Visualization of chromatin 
domains created by the gypsy insulator of Drosophila. 

/. Cell Biol., 162, 565-574. 



46. Muravyova,E., Golovnin,A., Gracheva,E., Parshikov.A., 
Belenkaya,T., Pirrotta,V. and Georgiev.P. (2001) Loss of 
insulator activity by paired Su(Hw) chromatin insulators. 
Science, 291, 495^198. 

47. Savitskaya,E., Melnikova,L., Kostuchenko,M., Kravchenko.E., 
Pomerantseva,E., Boikova,T., Chetverina,D., Parshikov.A., 
Zobacheva,P., Gracheva.E. et al. (2006) Study of long-distance 
functional interactions between Su(Hw) insulators that can 
regulate enhancer-promoter communication in Drosophila 
melanogaster. Mol. Cell. Biol., 26, 754-761. 

48. PhillipsJ.E. and Corces.V.G. (2009) CTCF: master weaver of the 
genome. Cell, 137, 1194-1211. 

49. LeeJ.S., Kim,Y., Kim,I.S., Kim,B., Choi,H.J., LeeJ.M., 
Shin.H.J., KimJ.H., KimJ.Y., Seo,S.B. et al. (2010) Negative 
regulation of hypoxic responses via induced Reptin methylation. 
Mol. Cell, 39, 71-85. 

50. Trinklein,N.D., Aldred.S.F., Hartman.S.J., Schroeder,D.I., 
Otillar.R.P. and Myers.R.M. (2004) An abundance of 
bidirectional promoters in the human genome. Genome Res., 14, 
62-66. 

51. Li,Y.Y., Yu,H., Guo,Z.M., Guo,T.Q., Tu,K. and Li,Y.X. (2006) 
Systematic analysis of head-to-head gene organization: 
evolutionary conservation and potential biological relevance. 
PLoS Comput. Bio/., 2, e74. 

52. Birney.E., Stamatoyannopoulos.J.A., Dutta.A., Guigo.R.. 
Gingeras.T.R., Margulies.E.H., Weng,Z., Snyder,M., 
Dermitzakis,E.T., Thurman.R.E. et al. (2007) Identification and 
analysis of functional elements in 1 % of the human genome by 
the ENCODE pilot project. Nature, 447, 799-816. 

53. Kent.W.J., Sugnet,C.W., Furey,T.S., Roskin,K.M., Pringle.T.H., 
Zahler,A.M. and Haussler,D. (2002) The human genome browser 
at UCSC. Genome Res., 12, 996-1006. 

54. Rosenbloom,K.R., Dreszer,T.R., Pheasant, M., Barber.G.P., 
Meyer,L.R., Pohl.A., Raney.B.J., Wang.T., Hinrichs.A.S., 
Zweig,A.S. et al. (2010) ENCODE whole-genome 

data in the UCSC Genome Browser. Nucleic Acids Res., 38, 
D620-D625. 



